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ABSTRACT 

We present a spectroscopic classification system for M-dwarfs and M-subdwarfs 
based on quantitative measures of TiO and CaH features in the region AA6200 — 7400A 
Our sample of cool stars covers the range from solar metallicity stars to the most 
extreme subdwarfs known. Using synthetic spectra computed by Allard and Hauschildt 
(1995), we derive metallicities for the stars. Stars are classified as dwarfs (M V), 
subdwarfs (sdM), or extreme subdwarfs (esdM). These classifications correspond to 
[m/H] ~ 0.0, -1.2, and -2.0 respectively. Our metallicity scale agrees with theoretical 
HR diagrams and HST globular cluster measurements. We discuss some nearby 
subdwarfs of particular interest in light of our metallicity scale. 



1. Introduction 

The vast majority of stars are M-dwarfs, main sequence stars whose spectra are dominated 
by molecular absorption. They have lifetimes much greater than the age of the universe which 
makes them an important fossil record of Galactic history. Their potential is largely unrealized 
because investigations of their properties have been hampered by the complex absorption spectra 
of diatomic and triatomic molecules. In particular, theoretical model atmospheres face serious 
difficulties; as a result, traditional methods of determining abundances from high resolution 
spectra of weak atomic lines are inapplicable, requiring the use of other techniques. 

A fundamental astronomical tool is a classification system that spans the range of observed 
properties, allowing an estimate of the effective temperature (T e jj), luminosity, and abundance 
{[m/H]) from a spectrum by comparison to standard stars with known properties. The system 
also should also allow identification of rare objects with unusual properties. The MKK dwarf 



spectral sequence flMorgan et al. 1943[) extended to M2 V and only included Population I objects. 



As progressively cooler stars have been discovered the classification system has been extended to 



M6.5 by Boeshaar (1976) and M9 flKirkpatrick et al.l99l[ , hereafter KHM; see also Bessell 1991). 



Observations were made partially at the 60-inch telescope at Palomar Mountain which is jointly owned by the 
California Institute of Technology and the Carnegie Institution of Washington 
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The KHM system uses features in the wavelength ranges AA6950 — 7500A and AA8400 — 8950^1 
observed at ~ 18A resolution or better, and Henry et al. (1994, hereafter HKS) have applied 
it to cool stars within 8 parsecs. A slightly different approach was taken by Reid et al. (1995, 
hereafter RHG) and Hawley et al. (1996, hereafter HGR) who used measurements of the 7100 
TiO bandhead at higher resolution (~ 3^4) to classify most of the known M-dwarfs within 25 
parsecs. The standards from KHM and HKS were used to place the RHG observations on the 
standard system. 

The spectral classification system is thus now well-defined for the near solar abundance M- 
dwarfs of the Galactic Disk. The situation for low metallicity stars, traditionally called subdwarfs, 
is much more confused. Although many spectra of individual stars have been published, there 
is no consistent M-subdwarf classification system. Examples of metal-poor M-subdwarfs were 
spectroscopically identified as long ago as Joy (1947), although "later investigators had difficulty in 
recognizing his criteria" ( port 19*65 ). Mould and McElroy (1978) discussed "old disk subdwarfs" 



which were less metal poor than other subdwarfs on the basis of TiO and CaH indices. Ake and 
Greenstein (1980), following a spectroscopic survey of high velocity stars, published spectra of 
four "extreme subdwarf M stars" which appeared to have "extreme metal deficiency" compared 
to the usually recognized M-subdwarfs. Similar stars were identified spectroscopically in a search 
for nearby white dwarfs ( Liebert et al. 1979| ), a search for Population II halo stars ( Hartwick et 



al. 1984 ), and a survey of cool M-dwarfs ( |Bessell 1982| ), all targeted at faint, high-proper motions 
stars. Recently, trigonometric parallaxes of 17 extreme subdwarfs have been measured, confirming 
their subluminosity QMonet et al. 19*92 , hereafter M92). Thus objects selected by different criteria 



can be called subdwarfs by different authors. M92 also found that there is a gap between their 
extreme subdwarf sequence and the "less extreme" subdwarf and disk sequence in the My vs. 
V-I HR diagram - their preliminary interpretation was that this gap represented a real lack of 
"intermediate" metallicity stars. In any case, the physical properties of the various subdwarf 
types remain in doubt; for example, M92 argue that their extreme subdwarfs have [m/H] ~ —1.7 
whereas Eggen (1996) suggests they have [m/H] ~ —2.5 to —3.5. 

Here we present a self-consistent set of spectroscopic observations of cool metal-poor stars 
drawn from a variety of sources. Selection of our subdwarf candidates and the data reduction 
procedures are discussed in Section ^. In Section ||], we discuss empirical molecular bandstrengths 
and suggest a two-dimensional classification system. In Section we compare our spectra to 
model atmospheres and deduce metallicities. In Section || we discuss the use of color-color and HR 
diagrams. Some notable individual stars are discussed in Section |(| The results are summarized 
in Section ^. 



2. Observations 
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2.1. Sample Selection 

Identification of a metal-poor or halo population can be difficult and ambiguous. It is well 
known from studies of hotter (F, G, and K) stars that there is no one-to-one correspondence 
between kinematics and metallicity because there is overlap between the properties of the 
disk and halo components (e.g., Mihalas and Binney 1981). Even among local stars with 
v tan > 100 km s _1 , the halo (Population II) is outnumbered approximately ten to one flSchmidt 



1975 ) by high velocity, slightly metal-poor disk stars (the Intermediate Population II, reviewed 



by [Majewski 1993 ). Despite the name Intermediate Population II (IPII), these stars have a mean 



metallicity ([m/H] ~ —0.6) close to that of the disk ([m/H] ~ 0.0), although the exact abundance 
distribution of the IPII is difficult to determine flCarney et al. 1989| ). In contrast, a conservatively 
selected field halo sample ( Laird et al. 19"88[ ) (restricted to stars that have vtan > 220 km s _1 



or V < -220 km s" 1 Q) has a distribution in [m/H] that peaks at -1.7 with FWHM ~ 1.2 dex, 
and includes tails that extend down to very low abundances (8% of the stars have [m/H] < —2.5) 
and up to IPII-like abundances (9% have [m/H] > —1.0). We therefore a priori expect that 
kinematically selected samples with loose selection criteria (e.g., vt an > 100 km s _1 ) will include 
many IPII stars that have [m/H] > — 1, and we also expect to see a large range of metallicities 
even among "true" halo stars. 

We have chosen our objects from a number of sources; however, virtually all have been 
identified in proper motion surveys and appear in the Luyten LHS catalog ( Luyten 197S| ) and 



most are in the the Lowell Proper Motion Survey flGiclas et al. 1971| ). These surveys give proper 
motions, photographic magnitudes, and low-precision color indices which are by themselves not 
adequate to isolate a halo or metal-poor sample. However, many late-type candidate subdwarfs 
have been identified in followup surveys. Since we use a variety of sources, selection criteria are 
ill-defined but typically depend on the star's velocity, color, or spectral features. We particularly 
favor stars with measured trigonometric parallaxes. These effects are not important for this work 
since we are not setting out to measure statistical quantities. Of particular note are the Schmidt 
(1975) complete sample of stars with fi > 1.295" yr _1 and m pg < 15.95, the Greenstein (1989) 
sample of cool halo stars, and the M92 and Ruiz and Anguita (1993, hereafter RA) samples of CCD 
parallax stars. Other objects have been chosen from lists of unusual stars in photometric followups 



QDawson and Forbes 1989] pawson and Forbes 1992] , [Reid 1982| , [Leggett 1992|) Finally, it should 
noted that the Giclas catalog only goes to m pg ~ 17 and is incomplete for m pg > 16. Schmidt 
(1975) showed that the apparent lack of M pg > 13.5 [My ~ 12) is due to this apparent magnitude 
limit. The Luyten searches reach mji ~ 20 but are incomplete both for mn > 18 and /i > 2.5"/yr 
( Dawson 1986| ). The latter limit implies that stars with typical halo tangential velocities of 220 



and 300 km s 1 would not have been detected within 18 and 25 parsecs respectively. 

In addition to the candidate subdwarfs selected above, we have utilized the RHG observations 



2 We use the standard notation of (U,V,W) for the space velocity components; note that U is positive towards the 
Galactic Center (I = 0,6 = 0). 
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of 1700 M-dwarfs from the preliminary Third Catalog of Nearby Stars ( pliese and Jahreiss 1991 



CNS3) as a reference solar abundance sample. We note particular use of two subsets from this 
sample. First, we use the stars within eight parsecs, selected to be either single or well separated 
from their companions, which outline the detailed structure in the HR diagram (e.g., Gizis and 
Reid 1996). Distances were adopted from RHG, but the stars were required to have accurate 



trigonometric parallaxes with Lutz-Kelker (1973) corrections less than 0.1 magnitudes. Second, 
we also use the nearby stars with vtan > 100 km s^ 1 . 



2.2. Spectroscopic Observations and Data Reduction 

Spectra were obtained at the Palomar 60 in. telescope, the Hale 200 in. telescope, and 
the Las Campanas Du Pont 100 in. telescope. We used the G-mode of the Palomar 60 in. 
spectrograph ( McCarthy 1985| ), a 1 arcsecond slit, and a 600 1/mm grating blazed at 6500A, 



yielding 1.5 A pix - . Candidate subdwarfs were observed in May 1994, September 1994, and 
January 1996 while some CNS3 stars were observed with the same setup in 1993 and 1994 (RHG). 
At the Hale 200 in. telescope we used the double spectrograph. In August 1995 with the blue 
camera was set to observe 6000 — 6900A and the red camera was set to 6700 — 8000A using 600 
1/mm gratings blazed at 4000A and IOOOOA respectively. In October 1995, a new red camera was 
installed in the double spectrograph, which we used in all subsequent runs to observe the region 
A6000 — 7400A at 1.4 A pix -1 . The 100 in. telescope observations used the modular spectrograph 
with a 1200 line grating blazed at 7500A The resolution at all telescopes was 3 — 4A. At all 
telescopes, neon and argon arcs were taken after each observation in order to eliminate the effects 
of instrument flexure. The data were extracted, sky subtracted, and wavelength calibrated using 
the FIGARO package. Flux standards ( |Gunn and Oke 1985 and Baldwin and Stone 1984 ) were 



used to set the data on an F u scale. Since observing conditions were usually non-photometric with 
seeing worse than 1 arcsecond, the fluxes are not absolute. 

Radial velocities were measured using the FIGARO cross correlation task SCROSS and 
corrected to a heliocentric frame using VHELIO. M-dwarf radial velocity standards were drawn 
from Marcy and Benitz (1986). For the Hale 200 in. telescope observations, the extreme subdwarf 
LHS 1174 was also used as a cross correlation standard. A radial velocity of —112 ± 2 km s _1 was 
adopted based upon a Keck HIRES spectrum (Reid, personal communication). Using either the 
extreme subdwarf or disk M-dwarfs did not make a significant change in the derived velocities, 
which have an uncertainty of ± ~ 20 km s~ . We have 14 stars in common with the more precise 
measurements of Dawson and De Robertis (1988, 1989) - we find a mean difference of 10 km s -1 
and a standard deviation of 26 km s~ . The three stars with differences of more than 30 km s -1 
are LHS 161, 205a, and 479. 
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3. Bandstrengths and a Subdwarf Spectroscopic Sequence 

Molecular features are typically very broad (tens or hundreds of Angstroms) and often are 
asymmetric. In M-dwarfs, there are few points that are relatively clear of absorption, but model 
atmospheres (Allard and Hauschildt 1995) show that even these are depressed from the "true" 
continuum. Since absolute measurements of the strength of features are therefore impossible, we 
instead measure molecular features with bandstrength indices defined as 

Rind — 

r cont 

where pseudo-continuum [F C ont) region and the feature wavelength limits are taken from 

RHG and listed here in Table [j]. In the case of the CaHl index, two sidebands are used to estimate 
the pseudo-continuum flux [S\ and 52). As in RHG, the spectral indices have an accuracy of 
±0.02 — 0.04. Note that our experience with these indices show that all except CaH2, which uses 
widely separated pseudo-continuum and feature regions, are insensitive to the flux calibration 
errors. CaH2 shows systematic offsets of up to 0.03 from observing run to observing run. The 
measured indices for our stars are given in Table |2[ 

The KHM M-dwarf spectral types provide an excellent shorthand description of the spectral 
properties of Galactic Disk stars - we find that over the range K7 - M6 the relations 

Sp = -9.64 x Ti05 + 7.76 (1) 

Sp = 7.91 x CaH2 2 - 20.63 x CaH2 + 10.71 (2) 

Sp = -18.00 x CaH3 + 15.80 (3) 

are accurate to ±0.5 subclass, where K5 and K7 are Sp = —2 and -1 respectively. As described 
by RHG, the CaHl index saturates at ~ M3, beyond which it does not show much temperature 
dependence. Although Ti05 was used by RHG to classify the stars, the two CaH relations are 
equally good spectral type indicators for near solar metallicity stars. 

Figure || shows the three CaH indices plotted against Ti05 for our candidate metal-poor 
stars and the eight parsec disk sample. Q Fifty stars lie significantly below the mean CaH vs 
Ti05 relations, a significant difference compared with only 5 of 1685 stars in RHG. There appears 
to be a fairly well defined sequence to the lower right of each figure - in Figure [j]b a high-order 
polynomial illustrates the cutoff adopted. We will call the stars below this cutoff the extreme 
M-subdwarfs. Note however that the region between the extreme M-subdwarf sequence and the 
disk M-dwarfs is well represented in our sample. 

The extreme subdwarf sequence, with the notable exception of LHS 453, has very weak but 
detectable TiO absorption at the resolution of our observations. However, even the coolest star, 



3 The relations defined by all of the M-dwarfs in the northern CNS3 are plotted in RHG's Figure 4. 
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LHS 1742a {My = 14.43), has only the TiO strength of a disk Ml star - the rest of the sequence is 
equivalent to K7 or K5 stars according to equation [[]. Indeed, LHS 1174 was classified "sdK7-Ml" 
by Liebert (1991, personal communication reported in M92). A TiO classification thus compresses 
the entire sequence of M-subdwarfs into only three subclasses even though Figure |l| shows there 
is as much or more variation in the CaH strengths as TiO variation in the disk M-dwarfs. Use of 
equations § or ||| gives spectral types of M5.5 (CaH2) or M8 (CaH3, extrapolated) for LHS 1742a. 

In many circumstances a simplified shorthand notation may be more useful than measurements 
of specific molecular features, even though Figure [l| shows that there is a full two-dimensional 
continuum in the properties of M-subdwarfs. Using the TiO and CaH indices as a guide, a 
sequence of low abundance M spectral standards can be defined. While this is the first subdwarf 
standard sequence to be defined, historical precedent supports using the prefix "sd" before the 



M spectral class. (Luminosity class VI should not be used for subdwarfs according to laschek 



and Jaschek 1987| .) In the past, LHS 64 has been classified sdMl (Joy 1947), very similar to the 



"sdK7-Ml" classification for LHS 1174. Unfortunately the use of "sd" alone throws away much 
information - although both types of subdwarfs above are significantly different from Population I 
stars, they also differ at least as much from each other, as seen in Table |2| and Figure [I]. Second, 
as discussed above, use of TiO alone would artificially suppress much of the variance among the 
extreme subdwarfs. 

We therefore propose the following classification system, which consists of both a prefix (M, 
sdM, esdM) and a numerical subclass. The classification system works as follows. First, determine 
if the star if a subdwarf of some kind using CaHl. Second, use CaH2 to determine if the star is a 
sdM or esdM. Finally, assign a numerical subclass. Each step is described in detail below. 

In the first step, stars that show significantly stronger CaH absorption for their Ti05 
strength compared with disk dwarfs are designed subdwarfs with the prefix "sdM" or "esdM". 
Quantitatively, we define significant as lying 0.07 below the mean Population I CaHl-Ti05 
relation of the RHG stars. If the other indices are unavailable, 0.06 below the RHG CaH2-Ti05 
or CaH3-Ti05 relations is equivalent. For the coolest stars (Ti05 < 0.49), CaH2 or CaH3 must 
be used since the disk CaHl relation is "saturated." The equations for these cutoffs are: 

CaHl < 0.695 x Ti05 3 - 0.818 x Ti05 2 + 0.413 x Ti05 + 0.651 (4) 

CaH2 < 0.968 x Ti05 3 - 1.358 x Ti05 2 + 1.315 x Ti05 - 0.033 (5) 

CaH3 < 0.639 x Ti05 3 - 1.199 x Ti05 2 + 1.161 x Ti05 + 0.307 (6) 

Stars that fulfill Equations |], [BL or || are either sdM or esdM. Equation ||] is plotted in 
Figure |]a. In the second step, the extreme subdwarfs are identified. A high order polynomial, 
shown in Figure [l|b, was used to interpolate between arbitrarily chosen points to define a 
cutoff in CaH2 vs. Ti05. The coordinates of the points, expressed as (Ti05,CaH2), are 
(0.0,-0.1), (0.433,0.234), (0.600,0.303), (0.800,0.456), (0.905,0.626), (0.977,0.788), and (1.020,0.914). 
The sdM lie above these points and the esdM lie below. The separation between the subdwarfs 
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and extreme subdwarfs is natural (at least for this sample) for the M stars, but at the hotter 
end the weakness of the features makes the sequences merge. At the cool end, the relation is 
not defined beyond Ti05 ~ 0.6, but the points given prevent the interpolating polynomial from 
inconveniently crossing the disk sequence. 

Finally, the numerical subclass is determined. The definition of the numerical subclasses is 
problematic. Larger numbers correspond to lower temperatures, but since effective temperatures 
are not well determined for cool stars it is impossible to ensure that sdM subclasses correspond to 
the same temperatures as M V stars or esdM stars. Some authors (e.g., Hamilton and Stauffci] 



1993J) have used pseudo continuum indices to assign types for late-type disk M-dwarfs, but we 
find that the pseudo continuum points available in our spectral range show more scatter than the 
bandstrength indices we are using. We adopt Equation ^ as the definition of numerical subclasses 
for the sdM and esdM because this feature shows the same range (0.2 <^ CaH2 1.0) for all 
metallicities. Use of the disk relation for the CaH3 feature (Equation ||) or the Ti05 feature 
(Equation |l]) leads to very large and very small ranges in the spectral type for the esdM. The 
subclass, however, is relatively uncertain when using only the CaH2 feature since the spectral 
type determination then depends upon only a single feature. We have therefore determined the 
following linear fits: 

Sp s dM = -16.02 x CaH3 + 13.78 (7) 

SpesdM = -13.47 x CaH3 + 11.50 (8) 

The numerical subclasses listed in Table |2] are the average of the CaH2 relation (Equation |2|) and 
the appropriate CaH3 relation (Equation ||[ 0, or ||). We continue to classify stars as "sd" or 
"esd" as early as K7. For earlier stars, we assign classes of only sdK (or sdG for stars with strong 
Ha absorption). The cooler sdK can be distinguished in Figure || where they have stronger CaHl 
relative to the K stars presented in RHG. These features are quite weak so the measurement errors 
are more important. The classifications of the K stars are therefore not as certain and should be 
regarded with caution - they are presented here to show the continuity with the M-subdwarfs. 
The new classifications are listed in Table @. Fi gures § and || show sequences of spectra for the 
subdwarfs and extreme subdwarfs respectively. 

The cool ends of the sequences are not yet well-defined, since there are as yet observations of 
only a few very cool spectroscopic subdwarfs. LHS 407 has been classified as sdM5.0. Although 
in Figure |l| it is well to the left of the esdM, it lies with the extreme subdwarfs in Figure |] 
and the HR diagram (Figure |7|). It is easily distinguished from LHS 1742a (esdM5.5), which 
is at present the coolest extreme subdwarf known. LHS 3061 is similar enough to LHS 1742a 
to also be esdM5.5, but it is slightly hotter according to its spectral indices and M92 V-I color 
(Ay-i = 0.13). LHS 3409, observed by RHG, is apparently more metal-rich than LHS 407, yet its 
indices are significantly different from the Population I, or even most stars with velocities typical 
of the IPII. We have therefore called it sdM4.5. LHS 3480 is a similar star. A considerably redder 
sdM is the star LHS 377 (shown by M92 to have My = 15.66). Its Ti05 index is equivalent to 
an M5.5 V star, but CaH2 gives the classification sdM7 (M92 report that Boeshaar and Liebert 
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call the star sdM5 in a private communication). Finally, we have also observed TVLM 832-42500 
which is subluminous in the I-K vs Mr diagram, which suggests it has [m/H] ~ —0.5 to —1.0 
( Tinney et al. 1995| ). It is at spectral type ~ M7 V where the indices break down; however, 
compared to the RHG stars of similar Ti05 strength it has slightly stronger CaH2 and CaH3, 
suggesting it is indeed slightly more metal poor, or at least spectroscopically different. Its Ti05 
and CaH2 absorption is slightly stronger than LHS 377 but the CaH3 absorption is much weaker. 
It thus appears to be more metal rich than LHS 377. A proper understanding of the distinction 
between cool sdM vs. esdM will require identification of more very cool stars from the LHS 
catalog. However, since we determine the M spectral subclass by the quantitative CaH2 index, the 
classification of M4.5 and M5.5 for the stars above is well-defined in our system. 

The relative strengths of bandstrength indices contains more information than just a simple 
classification. HGR show that there are systematic differences in the TiO indices of dMe and 
dM stars. In Figure ||, we plot the ratio CaH3/CaH2 as a function of CaH2 strength for all the 
RHG stars and for the sdM and esdM. The sdM tend to lie below the Population I sequence. 
The extreme M-subdwarfs lie systematically well below the sequence. The trend is that the first 
band of CaH (CaH3) is stronger relative to the entire bandhead (CaH2) - similar to the behavior 
seen in TiO for late type M-dwarfs. Since the deeper band must reflect conditions higher in 
the atmosphere, the difference might be caused by a slight surface gravity difference and/or a 
lower temperature in the CaH3 formation region compared to sdM and dM. Detailed modeling is 
necessary to understand this difference. 

4. Model Atmosphere Fitting 

Interpretation of the observed spectra in terms of the physical parameters effective 
temperature (T e ff) and metallicity ([m/H]) requires the use of model atmospheres. In principle, 
surface gravity is also a relevant parameter; however, both theory and observations indicate dwarf 
stars have logg ~ 5.0 ± 0.2 (Leggett et al. 1996) and we therefore restrict analysis to log g = 5.0 
models. We use the Extended Model grid computed by Allard and Hauschildt (1995, hereafter 
AH) which assumes local thermodynamic equilibrium (LTE). 

The theoretical uncertainties and limitations of the atmosphere models have been extensively 
discussed by AH. Although there are uncertainties in the appropriate treatment of the physics 
such as non-LTE effects, the treatment of molecular opacities (usually the Just Overlapping Line 
Approximation is used in these models), and the treatment of convection, the primary limitation 
is simply the lack of both laboratory data and theoretical computations of the electronic band 
oscillator strengths. In particular, no data exist for the CaH molecule which is extremely important 
in our spectral region. In Figure ||, we show a synthetic spectrum based upon a solar metallicity 
model and a representative disk star, the KHM M 4.0 V standard Gl 402. (The treatment of the 
TiO molecule in solar metallicity models has recently been improved by Allard and Hauschildt in 
the "NextGen" model). The spectra are in good general agreement, although the specific strengths 
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of the molecular features are incorrect, as expected. In particular, the CaH feature at ~ 6950^1 is 
too strong in all models compared to the observed spectra. However, AH show that their models 
capture the general overall behavior of metal-poor stars - with decreasing metallicity the hydride 
bands and atomic lines both increase in strength - the latter due to an increase in gas pressure. 

Figures [2] and |3| clearly show the expected increase in the atomic line strength with increasing 
hydride strength. Until accurate data or calculations are available for all the important molecules, 
but especially the hydride bands, any comparison between observations and theory will remain 
uncertain. Since such information is not likely to be available for many years, the uncertainties of 
this study will not reduced in the near future. Other authors have compared the AH models to 
cool solar abundance and metal-poor dwarfs using broad band colors and low-resolution optical 
and infrared spectra (Dahn et al. 1995; Lcggett et al. 1996| ) or higher resolution infrared spectra 



( Jones et al. 1996]) . Our study is most sensitive to variations in CaH and TiO which dominate our 



more limited wavelength region. It should also be noted that the metal poor models are computed 
assuming scaled solar abundances. If [O/C] is enhanced in the Population II M-subdwarfs, as it is 
in hotter Population II stars, it will affect the strength of TiO and H2O absorption ( [Mould 1976 ). 



Due to the uncertainties in the molecular data, the AH models cannot be used to predict the 
bandstrengths measured in Section |3], nor are they able to reproduce the CaH3/CaH2 relation in 
Figure ||. Instead, we compare the spectral region (AA6200 — 7300) by minimizing the least squares 
difference (x 2 ) between the observed spectrum and the models. We first shift the observations 
to rest velocity and vacuum wavelengths, convert to F\, and rebin to 2A steps. The models 
are convolved with a a = 2A Gaussian to approximate the instrumental resolution (in practice, 
this has negligible effects since the synthetic spectra often show 50A wide "features" due to 
the computation technique). The extended models have been computed for steps of 0.5 dex in 
metallicity between [m/H] = 0.0 and [m/H] = —4.0. We consider WOK steps for the temperature 
range 2500/T < T e ft < 4000-ftT (higher temperatures are not available and lower temperatures are 
irrelevant to the stars under consideration here). For temperatures that were not calculated (T e ff 
= 2600, 3100, 3400, 3600, 3800, and 3900 K), we approximate the synthetic spectra by linearly 
interpolating between the two nearest temperatures. 

The best fit model parameters for the subdwarf candidates are listed in Table ||. We have 
also fit the eight parsec sample of single disk M-dwarfs; below 3800 K, the spectra are best 
matched by solar metallicity models, as expected. Above 3800 K, the fits are poor because the 
molecular features are too weak to constrain the fits well. We therefore report only those fits with 
T e ff < 3800 K. The uncertainties from the fitting procedure are the grid spacing of ±0.5 in [m/H] 
and ±100 K in T e ff. The systematic uncertainties are unknown, but are dominated by the errors 
in the synthetic spectra discussed above. We note that different temperature scales for the disk 
M-dwarfs have historically had systematic differences of up to 300 K (e.g., |Leggett et al. 1996Q 



and our fits are subject to all the same uncertainties. Compared to Leggett et al. (1996), our fits 
give somewhat higher temperatures (AT e jf = 150K for LHS 377, AT e jj = 400 K for LHS 57). 
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The extreme subdwarfs have best fits at [m/H] ~ —2. Figure ^ shows one representative fit.. 
The coolest extreme subdwarf, LHS 1742a, has an effective temperature of 3300 K. A difference of 
one extreme subdwarf spectral subclass derived from CaH2 corresponds to an approximately 100 
K difference in temperature. Some authors (e.g., Eggen 1996) have argued that the the extreme 
subdwarfs are likely to have [m/H] = —2.5 to —3.0, whereas our least-squares fitting of model 
atmospheres gives metallicities a factor of 10 higher. The model atmospheres show that the most 
obvious qualitative difference is that ~ 7100^4 TiO features for the -2.5 and -3.0 models should be 
completely absent. The observed spectra, however, show that although the feature is weak it does 
exist in our extreme subdwarfs (except LHS 453), supporting the model atmosphere fits. LHS 
453 completely lacks the Ti05 feature (the TiO feature at 6200A is also missing) and therefore 
is more metal poor than the other esdM. The best fit is [m/H] = —2.5. For higher metallicities 
([m/H] = —1.0) the TiO feature remains strong which is qualitatively seen in the stars classified 
sdM. Indeed, the sdM have best [m/H] fits of -1.0 to -1.5. For the stars that were not classified 
subdwarfs, the model fits give [m/H] of 0.0 to -1.0. 

In summary, the spectroscopic classification system of Section |3] can be used to derive 
metallicities based upon our fits to AH synthetic spectra. Stars classified as extreme M-subdwarfs 
(esdM) have [m/H] of — 2.0=L —0.5, whereas sdM have a metailicity of — 1.2±0.3. Stars that do not 
show significant differences from the RHG disk sample in their CaH indices have [m/H] > —1.0. 
When model atmospheres are able to accurately reproduce the CaH and TiO bandheads more 
accurate metallicities should be possible. 



5. HR diagrams 

HR diagrams are an important tool for the study of cool stars. Observational HR diagrams 
are a critical test of stellar structure theory, provided that the transformation between theoretical 
temperatures and luminosities and observed quantities (colors and absolute magnitudes) can be 
made. In practice there has been disagreement over the Population I temperature scale ( [Leggett 



et al. 1996| ) and relatively little discussion of the temperature of Population II subdwarfs. The HR 



diagram can also be used to predict absolute magnitude from observed color and thus is used in 
interpretation of star count data and derivation of the luminosity function. Finally, derivation of 
the mass function from the luminosity function of metal-poor stars depends upon the quality of the 
theoretical mass-luminosity relations and evolutionary tracks since there are no cool Population II 
binaries with mass determinations. 

Table |] lists optical photometry and trigonometric parallaxes from the literature for the stars 
in our sample. We have adopted the Cousins R and I system and where necessary have used the 
color transformations compiled by Leggett (1992). Because the photometry is from heterogeneous 
sources of varying quality, in some cases the colors may be incorrect by as much as 0.1 magnitudes 
but the typical accuracy should be about 5% ( Leggett 1992| ). In addition, the transformations 



have been derived for near solar metailicity stars and may be inaccurate for subdwarfs. Absolute 
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parallaxes were taken from the preliminary Yale Parallax Catalog ( |van Altena et a!.1991| ) when 
available. In the case of the CCD relative parallaxes reported by RA, we added a correction of 0.5 
milliarcseconds to the parallaxes to transform them to absolute measurements. This approximately 
matches the correction applied by M92 for similar magnitude stars and the uncertainty in this 
correction is much less than the uncertainties in the relative parallaxes. Although the measured 
parallax is the best estimator of the distance to any particular star, there is a statistical bias to 
underestimate the distance and hence underestimate the the mean luminosity of the sample stars. 
Lutz and Kelker (1973) derived statistical corrections to the observed absolute magnitudes as a 
function of a/ir and showed that only parallaxes with a/ir < 0.2 are useful. However, selecting the 
stars by proper motion reduces the statistical correction ( Hansen 19781 ) . Since the selection of our 
sample is not well-defined, we do not apply any statistical correction to the absolute magnitudes, 
but we list the Lutz-Kelker corrections (derived from the approximation given by [Hansen l979| and 
defined as M true = M b s + A-lk) as a guide to likely biases in the sample. The V-I, My diagram is 
shown in Figure 0. The eight parsec single stars are plotted with Lutz-Kelker corrections applied, 
but the sample is restricted to stars with corrections less than 0.1 magnitude. We also plot the 
B-V version of this diagram in Figure H including the bluer K sub dwarfs in the sample and the 
disk G and K dwarfs within eight parsecs from the CNS3. It is well known that B-V colors are 
poorly suited for estimating absolute magnitudes for the disk M-dwarfs (in the disk sequence My 
changes by 4 magnitudes in the range 1.4 < B — V < 1.6.) In our diagram, the sdM and esdM 
sequences cross the disk sequence at B-V ~ 1.5, and the reddest sdM and esdM actually lie above 
the disk sequence for B-V > 1.65. 

Optical and near-infrared color-color diagrams for M-dwarfs have recently been extensively 
discussed by Leggett (1992). In Figure || we show the V-I vs. B-V diagram for the objects we have 
classified. The well-known tendency for metal-poor M-subdwarfs to have redder B-V at a given 
V-I (e.g., Mould and McElroy 1978| ; pahn et al. 1995| ) is clearly evident, for V - I > 1.7. The 



offset appears to be strongly related to metallicity, as the esdM are approximately 0.3 magnitudes 
redder than the disk stars in B-V with our sdM lying in between the two sequences. A V-I vs. 
B-V diagram is therefore useful in identifying red subdwarfs and extreme subdwarfs. However, 
giants also lie redwards in B-V at a given V-I ( Weistrop 1977|) . Indeed, Reid (1982) suggested on 



the basis of optical (BVRI) colors that Sm 183 might be a subdwarf similar to Kapteyn's star, 
but also noted the possibility that it could be a giant. In our spectrum, the CaH indices are weak 
compared to TiO, implying that it is a low surface gravity giant. 



5.1. Tests of the Spectroscopic Metallicity Scale 

Theoretical structure predictions and HST globular cluster observations provide checks on 
our spectroscopic metallicity scale. We compare with broadband colors and magnitudes for 
[m/H] = 0.0, -0.5, and -1.5 computed with stellar interior models that use the Allard and 
Hauschildt model atmospheres as outer boundary conditions ( [Baraffe et al. 1995|) . The tracks are 
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shown in Figure |T^. Baraffe et al. (1996) did not compute [m/H] = —1.0 models, but they did 
note that for [m/H] = —2.0 the models are shifted blueward by ~ 0.2 magnitudes with respect 
to the -1.5 track. Also shown is an even more recent set of solar metallicity models computed 
with the latest Allard and Hauschildt NextGen models (Baraffe and Chabrier 1996). The results 



dramatically illustrate the importance of boundary conditions - at 0.3Mq (My = 11.31, V-I=2.19, 
T e ff = 3403), the newer model is 201 K hotter, 1.21 magnitudes brighter in My, and 0.63 
magnitudes bluer in V-I. As a result, the newer models fail to match observations (suggesting that 
they would be a poor choice to use for the My-Mass relation). In our HR diagram, the earlier 
set of computations give better agreement with the Population I sequence, although none of them 
convincingly duplicate the structure seen. As it stands, these models support the identification 
of the esdM with [m/H] ~ —1.5 to —2.0 and the sdM as more metal poor than -0.5. Given 
the change of the theoretical models in less than one year, effective temperatures, metallicities, 
and masses derived from evolutionary models are clearly very uncertain - but the use of Allard 
and Hauschildt's extended model grid leads to consistent results for both spectroscopic and HR 
diagram metallicities. 

HST photometry of a stellar cluster yields both a luminosity function and a color-magnitude 
diagram. Santiago et al. (1996, hereafter SEG) provide piecewise linear fits to the HR diagram of 
the clusters M15 ([m/H] = -2.26), u Cen (-1.6), 47 Tuc (-0.6), NGC 2420 (-0.45), and NGC 2477 
(0.0), which are shown in Figure |ll| superimposed on our data (including the earlier subdwarfs). 
The globular cluster data reach only the brightest of the sdM and esdM. The offsets however 
appear to be consistent with the esdM having abundances between -1.6 and -2.26 and the sdM 
between -1.6 and -0.6, as deduced for our spectroscopic analysis. SEG note that the slopes of their 
CMD are steeper than the nearby field parallax subdwarfs and attribute this to either calibration 
errors in the HST photometry or real astrophysical differences. The considerable difficulties in 
transforming HST magnitudes to the Cousins system are discussed by Harris et al. (1991) and 
Holtzman et al. (1995). In particular, the HST F606W filter includes a substantial fraction of the 



Rc passband. However, the differences in slope of the main sequence are less evident in Figure |11 
than in the fits reported by SEG in the region of overlap. It appears likely that the confusion in 
the definition of subdwarfs led the the inclusion of progressively less metal-poor stars at redder 
colors in the Richer and Fahlman (1992) main sequence, in particular stars with V-I ^ 2.4 which 
we argue below are Intermediate Population II ([m/H] ~ —0.6) and which, in any case, are clearly 
much fainter and redder than the stars used to define the globular cluster slopes in the HST study. 
The fit to the NGC 2477 lower main sequence appears to lie above the local stars - perhaps due 
to color terms but perhaps also because of the inclusion of binaries or non-cluster members in the 
fits. 

Most of the cool extreme subdwarfs in this paper come from the M92 USNO CCD parallaxes. 
Our spectroscopic metallicity scale is in agreement with the M92 estimate that the mean [m/H] 
of extreme subdwarfs is ~ —1.7. Some authors have disputed the argument that the extreme 
subdwarfs represent the 'typical' halo metallicity of ~ —1.7. Richer and Fahlman (1992) argue 
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that the extremely hot kinematics (15 of the 17 M92 extreme subdwarfs have vtan > 275 km s _1 
and the 5 stars with known total space velocities have < vtot >= 427 km s _1 ) may point to 
their being more metal poor than the typical Population II stars. However, since there is no 
metallicity- velocity relation within the higher mass stars of the Population II halo ( [Norris 1986 ) 
this argument seems unconvincing. Note however that the Richer and Fahlman color-magnitude 
relation matches our sdM at V-I ~ 1.5 but matches the disk sequence by V — I ~ 3, which may be 
appropriate for their sample since their "Population II" luminosity function is in fact dominated 
by the IPII at the faint end ( Reid et al. 1996| ), Lacking a direct determination of metallicity, 
Eggen (1996) estimated values of -2.5 to -3.5 for the M92 extreme subdwarfs but considered these 
values "guesses only." 



5.2. Where are the cool sdM? 

M92 note that their HR diagram (their Figure 10) shows a remarkable gap between the 
extreme subdwarfs and the disk sequence. In their HR diagram, only one subdwarf appears 
between the two sequences (at about V-I of 2.2). In our sample, there is clearly a continuous 
distribution for stars with V-I< 2.2, leaving the gap for 2.2 < V-I< 2.8. There are three major 
solutions to the gap, each of which was addressed by M92. First, the gap may be an illusion 
caused by stellar astrophysics - a continuous distribution in metallicity may not necessarily 
correspond to a simple distribution in a given observational HR diagram. Second, the gap may be 
due to selection effects such that stars in the gap are not chosen for parallax programs. Finally, 
it may represent a real lack of stars with metallicities between the disk sequence and the extreme 
subdwarf sequence. We will reconsider each of these options. It is certainly true that bandpass 
effects and the properties of cool atmospheres are important in interpreting HR diagrams. In 
Figure ||, we have seen that the behavior of the B-V HR diagram is opposite to that of the V-I 
diagram - at the reddest colors, the subdwarfs lie above the disk sequence, and the esdM lie 
above the sdM. A second example of the importance of bandpass effects is in the V-I HR diagram, 
where it appears that for the very coolest disk stars (LHS 2924 and 2065) V-I turns around, 
and becomes bluer with decreasing temperature ( Bessell 1991| ; M92). Nevertheless, there is no 



evidence at present that the gap is caused by bandpass effects because M92 report the Allard 
model atmospheres do not produce such an effect. M92 consider that it is "very unlikely" that the 
gap is due to the USNO selection criteria, and therefore offer the preliminary interpretation that 
stars with intermediate metallicities are rare in the solar neighborhood. Since this result is quite 
surprising, we will next reconsider the selection effects that may entered this sample. 

There necessarily were important selection effects in the catalog, since as discussed in 
Section |||, even in a proper motion survey the majority of stars will not be extremely metal poor 
or halo stars. It is difficult to know exactly what the selection effects entering the catalog were, 
but early studies that emphasized the targets surely were influential. Some of the M92 stars were 
never studied before (according to SIMBAD), but others had already been noted as interesting. 
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Two of their stars (LHS 3382 and LHS 489) were pointed out by Ake and Greenstein (1980) as 
being spectroscopically much more metal poor than other known subdwarfs. Liebert et al. (1979) 
reported "sdM" spectra for five others (LHS 192, 197, 205a, 207, and 453) from a search for nearby 
degenerate stars - their Figure 2 presents a noisy spectrum of LHS 205a with strong hydride 
bands and a lack of TiO. Four others (as well as LHS 3382 again) were chosen by Hartwick et al. 
(1984) for having extremely large reduced proper motions, defined as 

H = M + 5 log iw/4.74 = m + 5 log p + 5 

The lack of the sdM is obvious in their reduced proper motion diagram (their Figure 5). In 
contrast, there is no obvious gap in the reduced proper motion diagram (based on B and R 
magnitudes from the POSS plates) of the entire LHS catalog presented by Dawson (1986, 
his Figure 7). One explanation of the very large tangential velocities of the M92 sample 
(< vtan >= 380 km s^ 1 ) is that there was a implicit selection effect favoring large H. If so, a group 
of subdwarfs ~ 1 magnitude more luminous but the same reduced proper motions would have 
< vtan >= 600 km s _1 . Since this is greater than the escape velocity of the galaxy (500 km s _1 , 
Carney et al. 1988| )), these subdwarfs could not enter the sample. Similarly, the (unexplained) 
lower limit of vt an = 275 km s _1 would correspond to a lower limit of vtan = 435 km s _1 . We 
therefore suggest that such a selection effect (as well as spectral selection of very weak TiO lines) 
would explain the lack of intermediate metallicity stars. As noted above, these stars are too faint 
to have been included in the more completely studied Giclas catalog which includes many bright 
early sdM. In any case, there are some likely LHS Catalog red sdM in the study of Hartwick et al. 
(1984) which do not yet have parallaxes - they identified both Class H stars (e.g., LHS 192, 1970, 
3259, 3382, and 3548), which we call esdM, and Class I stars (e.g., LHS 29 and 64) which we can 
identify as sdM stars. They found three faint class I stars: LHS 2533, 2630, and 3189. We confirm 
the latter as an sdM5.0 star. Further observations of these three stars would greatly increase the 
information available on cool sdM stars. 



5.3. Kinematics 

Using the trigonometric parallaxes in Table ||, the radial velocities in Table ^, and the 
LHS catalog proper motion values, we can compute the U,V,W components of galactic velocity 
(Table ||). The errors in radial velocity, ± ~ 20 km s _1 , are comparable to the uncertainties 
due to the trigonometric parallax errors. The kinematics of the sample are not representative 
of the Population II due to the uncertain selection criteria, but they can provide some 
information on the selection criteria used. Mean values and standard deviations (velocity 
dispersions) for the spectroscopic classes are listed in Table ^. The sdM, with the mean velocity 
component < V >= —202 km s^ 1 , and the esdM, with a mean galactic velocity component 
< V >= —287 km s _1 , have clear Population II kinematics (the mean V velocity for the halo is 
disputed but lies between -270 and —180 km s _1 according to [Majewski 1993| ). Our candidate 
halo stars that were classified as M V have both smaller < V >= —134 km s _1 and smaller U 
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and V velocity dispersions - more typical of the IPII and certainly much greater than the disk 
(< V >= -22 km s _1 ). These non-subdwarfs in fact have bandstrength indices similar to the 
majority of RHG nearby stars with Vt an > 100 km s _1 , which should be an IPII dominated 
sample. For both sdM and esdM, ajj and aw are somewhat greater and ay is somewhat less than 
the halo values of Norris (1986) as expected for a sample selected by large velocity but excluding 
near solar V values. Also given in Table ^ are the velocity dispersions for late M-dwarf RHG 
stars with TiOh < 0.4 and either weak CaHl (CaHl > 0.74) or strong CaHl (CaHl < 0.74), as 
well as the RHG values for nearby disk M-dwarfs. The strong CaHl stars have larger velocity 
dispersions, suggesting that they are older and supporting the idea they are more metal poor - 
but their velocities are much less than a Population II group. Thus although CaHl is useless as a 
temperature indicator for stars cooler than M3, it may prove to be useful metallicity indicator for 
late-type (M4 to M6.5) disk stars, particularly if CaH features can be reproduced by future model 
atmospheres. 

We have also computed total galactocentric velocity for the stars after correcting for the Sun's 
motion with respect to the Local Standard of Rest ( Mihalas and Binney 1981] ) and assuming a 



LSR circular velocity of 220 km s . All of the stars are easily bound to Galaxy - the fastest sdM 
(LHS 467, v Gal = 342 km s" 1 ) and the fastest esdM (LHS 3548, v Gai = 364 km s -1 ; LHS 453, 
v Gai = 451 km s _1 ) are comfortably less than the galactic escape speed of ~ 500 km s _1 ( Carney 
et al. 1988|). 



6. Notes on Individual Stars 

6.1. CM Dra (LHS 421) 

CM Dra (Gl 630.1) is a short period eclipsing binary with component masses 
M A = 0.2307 ± 0.0010 and M B = 0.2136 ± 0.0010M o and helium abundance Y A = 0.32 
and Yb = 0.31 ± 0.04 ( [Metcalfe et al. 1996 ). Based on the large W velocity and the observed 



low flaring rate and lacking a spectroscopic metallicity, [Metcalfe et al. 1996 interpret the system 



as Population II and therefore interpret the high helium abundance as possible support for a 
higher primordial helium abundance. Using the observed spectrum (taken from RHG), we find 
that the system not very metal poor, with a best fit model of [m/H] = 0.0, and does not show 
any CaH excess compared to Population I stars; however, infrared photometry ( Leggett 1992| ) 



suggests it is more metal-poor than most disk stars. We therefore conclude CM Dra formed from 
significantly enriched material. The CNS3 W velocity of -34 (-27 when corrected for the solar 
motion given in Mihalas and Binney 198l| ) implies that its maximum height above the plane is 



only ~ 350 pc ([Kuijken and Gilmore 1989| ), whereas starcounts show the true Population II does 
not become dominant until beyond 5 kpc ( Majewski 1993; ), so CM Dra can probably be considered 



an Intermediate Population II, or even Population I, star. 
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6.2. Barnard's Star (LHS 57) and Gl 299 (LHS 35) 

M92 have discussed the well known Barnard's Star (Gl 699) extensively. It lies perhaps 1.6 
magnitudes below their mean disk main sequence, yet it is well known that spectroscopically it 
is an M-dwarf, a result that this study confirms. We believe that the distance below the main 
sequence is "artificially" enhanced by the presence of the kink in the main sequence ( pizis and 



Reid 1996 Figure 5) which is not fit well by a line. At spectral type M4.5 the disk sequence, 
when restricted to only single stars with high quality parallaxes and photometry, shows a sharp 
1 magnitude drop. Hence Barnard's Star is only about 0.6 magnitudes below the disk sequence. 
It seems likely that Barnard's Star has a metallicity between —0.5 and —1.0. It is a weak outlier 
in TiO and CaH. The well known star Gl 299 (LHS 35) is quite similar in both its indices and 
absolute magnitude. We identify both as Intermediate Population II stars and as cool analogues 
to the spectroscopically unremarkable stars (e.g., LHS 301 and 376). The M4.5 V star LHS 3684 
which has an M92 parallax is an only slightly cooler counterpart of Barnard's Star in both the 
spectroscopic indices and the HR diagram. 



6.3. The sdMe (LHS 482 and 2497) and Unresolved Binaries 

Two subdwarfs in this study show noticeable Ha emission. Gl 781 (LHS 482) shows large 
radial velocity variations typical of a short period binary ( Joy 194"^ ). We have obtained echelle 



spectra that confirm the radial velocity variations but do not show lines from the secondary ( Gizis 



et al.1996 ). We therefore have not adjusted the absolute magnitude in Tabldl] for the secondary 



component. Gl 455 (LHS 2497) shows weak H alpha emission ( Hawley et al. 1996j ) and our 



recent echelle spectra have shown the system to be a short period double-lined spectroscopic 
binary flGizis et al.l996j ). Since the system is near-equal luminosity we have not adjusted the 
color but have subtracted 0.75 from the absolute magnitude in Table ||. The system is a strong 
outlier in Figure |], yet even with the adjustment it has an absolute magnitude and color similar to 
the subluminous IPII dM stars. The YPC parallax is based upon two determinations with good 
agreement and has an uncertainty of less than 10%. If it is not in error, perhaps the chromospheric 
activity affects the spectrum or magnitude of this interesting system. Unfortunately we have no 
other examples of late sdM that is comparable to this system. Young et al. (1987) argue that 
M-dwarf binaries with short periods (P <^ 5 days) must show emission. The observed radial 
velocity variations are consistent with this limit. Mass determinations would provide a important 
constraint on theories of metal poor low-mass stars, but the short periods imply small separations 
(a ^ 4 milliarcseconds) for these systems. No other sdM or esdM show Ha emission in this study, 
and indeed no extreme M-subdwarfs with emission are yet known. 

The sdMe discussed above are the only stars in the sample known to have unresolved 
companions. LHS 169, also known as Gl 129 but no longer within the nearby star catalog (CNS3) 
due to an improved parallax, has been reported to be spectroscopic binary. Three measurements 



-17- 



of the radial velocity (Joy 1947) found radial velocities of -139.6, -91.0 and —75.9 km s" 1 with 
probable error 2.1 km s _1 . The radial velocity variations have evidently never been confirmed or 
refuted since then, and Dawson and De Robertis (1988) suggest that the Joy (1947) velocities 
may be incorrect for a number of subdwarfs. LHS 64 has been reported to have velocity varying 
between -292 and -242 by Joy (1947) but we have obtained 3 echelle spectra which are consistent 
with no variations at all. Lacking any further data, we do not apply any correction to the absolute 
magnitude of either star for any companion. Presumably at least a few of the stars in our sample 
have unresolved companions that are still unknown, which could affect the magnitudes or radial 
velocities in a few cases. 



6.4. The most extreme subdwarf: LHS 453 



One important use of a classification system is identifying unusual stars. Our indices have 
allowed us to identify LHS 453 as even more metal poor than the other extreme subdwarfs. We 
would therefore expect this star to be below the other stars in the HR diagram, but the M92 
observations (V-I = 2.22, My = 13.08) place the star in the midst of the extreme subdwarf 
sequence. However, the reported mean error in My is still fairly large at ±0.20 magnitudes. 
Perhaps coincidentally, this most metal poor star also has the largest galactocentric velocity in the 
sample, though Dawson and De Robertis (1988) showed that despite the large tangential velocity 
(M92 give 466.9 ± 42.8) the star is bound to the galaxy. If the parallax is indeed reduced, the 
tangential velocity would also be decreased. Alternatively, if LHS 453 is a binary, its My would 
be fainter by up to 0.75 magnitudes. The best fit synthetic spectrum has [m/H] = —2.5, but we 
note that the synthetic [m/H] = —3.0 spectrum is almost identical. Although TiO is absent, LHS 
453 is not likely to be more metal poor than this because the observed atomic lines are prominent 
and similar in strength to those of LHS 205a (which is only slightly later type in V-I color and 
CaH indices). According to the models, the atomic lines should be much weaker or absent for 
[m/H] = —3.5 or -4 at all temperatures. 



7. Conclusions 

This paper presents spectroscopy of 50 late-type stars that are significantly more metal-poor 
than the Population I. Bandstrength indices are used to measure TiO and CaH features and show 
that many different metallicities can be distinguished in the Population II M stars, as expected 
from the metallicity spread seen in G-subdwarfs. The TiO feature at ~ 7100A is shown to put 
all the extreme subdwarfs into only a few spectral subclasses (K7 and M0). Instead, a CaH 
index which reproduces the Population I spectral classification system of Kirkpatrick, Henry, and 
McCarthy (1991) is used to classify stars. Based on the differences between the (weak) TiO and 
the CaH the metal poor M stars are divided into the simple categories of M-subdwarfs (sdM) and 
extreme M-subdwarfs (esdM) and spectral standards are presented (Figures || and|||). 
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Synthetic spectra computed by Allard and Hauschildt (1995) are fit to the observed spectra. 
Metallicities of [m/H] ~ —2 are derived for the esdM and ~ —1.2 for the sdM. These estimates are 
in agreement with the estimates of G89 and M92 but not those of Richer and Fahlman (1992) or 
Eggen (1996). HR diagrams based on literature observations are presented for dM, sdM and esdM. 
They are at least qualitatively in agreement with theoretical computations and HST observations 
of globular clusters. However, the former are sensitive to the boundary conditions used and the 
latter are likely to suffer from large color terms. The agreement of all three methods seems to 
indicate the metallicity scale is accurate to ~ 0.5 dex. 

Given this metallicity scale, we discuss a number of interesting subdwarfs and stars that have 
been called subdwarfs. The eclipsing binary system CM Dra is shown to be significantly enriched 
such that its composition is not likely to reflect the primordial helium abundance. The "old disk 
subdwarfs" such as Gl 299 and Gl 699 are probably not as subluminous as is usually thought - 
taking into account the detailed structure of the HR diagram, rather than simply fitting a straight 
line, they are only ~ 0.6 magnitudes subluminous rather than the usually cited ~ 1.2 magnitude. 
This naturally explains the only slight spectral differences, which are shown to be similar to 
differences of hotter M-dwarfs that are also ~ 0.6 magnitudes subluminous. Finally, we show that 
LHS 453 is more metal poor than the other extreme subdwarfs with [m/H] ~ —2.5 to -3. 

I would like to thank Neill Reid for helpful comments and help with the observations. The 
referee, Todd Henry, provided comments which substantially improved this paper. I am indebted 
to France Allard and Peter Hauschildt for making their synthetic spectra available. I am grateful 
for Greenstein and Kinsley Fellowships and partial support through NASA grant GO-5353.0-93A. 
This research has made use of the Simbad database, operated at CDS, Strasbourg, France. 



Table 1. Spectroscopic Indices 



Band SI W S2 



TiO 5 7042-7046 7126-7135 

CaH 1 6345-6355 6380-6390 6410-6420 

CaH 2 7042-7046 6814-6846 

CaH 3 7042-7046 6960-6990 



Table 2. Spectroscopic Observations 



LHS Name Ti05 CaHl CaH2 CaH3 v rad class Source 



12 


G 003-036 


0. 


.894 


0. 


.852 





,746 


0. 


.878 


25. 


.1 


sdMO.O 


60 


20 


GJ 1062 





.677 


0. 


.708 





,499 


0. 


710 


-104. 


.3 


sdM2.5 


RHG 


29 


Kapteyn's 





,810 


0. 


.778 


0. 


,594 


0. 


.792 


242. 


.8 


sdMl.O 


100 


42 


Ross 451 


0. 


.925 


0. 


.871 


0. 


.757 


0. 


.885 


-148. 


,0 


sdMO.O 


60 


55 


GJ 1200 





.439 


0. 


.816 





.420 


0. 


.705 


-21 


.2 


M3.5 V 


60 


57 


Barnard's 





.394 


0. 


.752 





.394 


0. 


.652 


-120. 


.2 


M4.0 V 


60 


61 


Gl 817 





.645 


0. 


.823 


0. 


.582 


0. 


.799 


25. 


.9 


M1.5 V 


60 


64 


Wo 9722 





.773 


0. 


.731 


0. 


.572 


0. 


.754 


-235. 


.5 


sdM1.5 


RHG 


104 


G 030-048 


1 


,009 


0. 


911 


0. 


.854 


0. 


.938 


-173. 


.8 


esdK7 


60 


156 


G 004-029 


0. 


.655 


0. 


.648 


0. 


.445 


0. 


.659 


-42. 


.5 


sdM3.0 


100 


161 


G 075-047 


0. 


,964 


0. 


739 


0. 


.515 


0. 


.751 


-35 


,1 


esdM2.0 


100 


169 


G 005-022 


1 


.008 


0. 


.879 


0. 


.862 


0. 


.951 


-119. 


.2 


esdK7 


60 


170 


G 078-026 


1 


,007 


0. 


.986 





.931 


0. 


.961 


-183. 


.2 


sdK 


60 


173 


G 038-001 


0. 


.961 


0. 


.942 





.868 


0. 


938 


-159. 


.2 


sdK7 


60 


174 


G 037-040 





,877 


0. 


.817 





.686 


0. 


.833 


-225. 


.2 


sdM0.5 


60 


178 


G 079-059 


0. 


.674 


0. 


.766 





.551 


0. 


.760 


-59. 


.7 


sdM1.5 


60 


182 


G 095-059 





.975 


0. 


.842 


0. 


.735 


0. 


.843 


-238. 


.1 


esdMO.O 


60 


185 


G 007-017 


0. 


.957 


0. 


737 


0. 


.659 


0. 


.795 


37. 


,7 


esdM0.5 


60 


192 


LP 302-31 





.946 


0. 


.710 





.603 


0. 


.768 


79. 


.9 


esdMl.O 


200 


205a 


LP 417-44 





,889 


0. 


.512 


0. 


.347 


0. 


513 


-76. 


.4 


esdM4.5 


200 


211 


G 099-033 


0. 


.796 


0. 


.797 


0. 


.703 


0. 


.837 


-129. 


.2 


sdMO.O 


100 


216 


G 105-023 





.725 


0. 


.743 





.496 


0. 


732 


392 


.5 


sdM2.0 


60 


218 


G 103-046 





.500 


0. 


.759 


0. 


.413 


0. 


.686 


-20. 


.4 


M3.5 V 


100 


236 


G 251-044 





.980 


0. 


.944 


0. 


.879 


0. 


931 


66 


.8 


sdK7 


60 


254 


LP 666-11 





.178 


0. 


.678 


0. 


.230 


0. 


517 


20. 


.2 


M6.5 V 


100 


272 


LP 788-27 





.678 


0. 


.676 





.441 


0. 


.669 


264. 


.4 


sdM3.0 


100 


276 


G 117-061 


0. 


.686 


0. 


.888 





623 


0. 


.828 


46. 


.9 


M1.0 V 


60 


301 


GJ 1146 





.481 


0. 


.773 


0. 


.431 


0. 


.683 


94. 


.4 


M3.5 V 


60 


307 


G 176-040 


0. 


.910 


0. 


763 


0. 


.671 


0. 


.837 


-64. 


.8 


sdM0.5 


60 


320 


G 011-035 





.620 


0. 


.756 


0. 


.512 


0. 


.751 


81, 


.1 


sdM2.0 


60 


343 


G 061-021 


1 


.001 


0. 


953 





.901 


0. 


.937 


165 


.5 


sdK 


60 


364 


G 165-047 





.979 


0. 


699 


0, 


.604 


0. 


.740 


12. 


,8 


esdM1.5 


60 


375 


LP 857-48 





,829 


0. 


.499 


0. 


.372 


0. 


.547 


181. 


.2 


esdM4.0 


60 


376 


G 135-067 





.448 


0. 


.795 





.423 


0. 


.682 


-19. 


.6 


M3.5 V 


60 


377 


LP 440-52 





.232 


0. 


.567 


0. 


.205 


0. 


396 


179. 


,7 


sdM7.0 


200 


407 


LP 803-27 





.601 


0. 


.584 


0. 


,351 


0. 


.528 


-173. 


.1 


sdM5.0 


60 
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Table 2 — Continued 



LHS 


Name 


Ti05 


CaHl 


CaH2 


CaH3 


Vrad 




class 


Source 


41 n 

~t 1 VJ 


n nin-Di8 

VI VJXU VJXO 


794 

\J. 1 ill 


D 899 


D 61 9 


n 814 


i n 

XVJ. 




mi n v 

1VXX.VJ V 


6fl 

VJVJ 


418 


G 1 38-f)25 

V.I lOU VJZj>J 


941 

VJ . JtI 


D 9D8 

VJ . sJVJO 


n 814 

V J . O X 'X 


n 9i n 

\J . -J XVJ 


-14. 


4 


K7 V 

XV 1 V 


6fl 

VJVJ 


425 


O 1 38-D59 

V.I iOU VJOiJ 


f) 646 

VJ . VJtVJ 


762 

VJ . 1 VJ Zj 


D 5D3 


D 738 


-49 


< 

O 


srlM2 D 

oVXIVXZj . VJ 


6fl 

VJVJ 


453 

x 'J (J) 


LP 139-14 


1 070 

X . VJ t VJ 


514 


407 


590 


32 

O^ . 


1 


esdM3 5 

CoVXIVXO. KJ 


200 


460 


GT 1225 

V.I t J x ^ ^ o 


337 

VJ. 1 


736 

VJ . 1 OVJ 


321 


603 


-53 


2 


M5 V 

1VXO.VJ V 


60 

VJVJ 


467 

i VJ i 


O 091-093 

V.I Ua*1 IJZjO 


D 966 

VJ. JJUVJ 


fl 898 

VJ .0i70 


D 786 


D 897 


1 83 

xoo . 


VJ 




6D 

VJVJ 


479 


d 149-D59 

V.I ItIi ua^ 


D 799 

VJ. 1 


fl 768 

U. 1 VJO 


n 6i n 


n 8D1 

U.OU X 


-75 

( o. 


4 


oVllVXX . VJ 


6fl 

VJVJ 


482 


Gl 781 

V.I 1 | O 1 


774 


n 718 


D 558 


D 746 


-1 D7 

XVJ 1 . 


fi 

VJ 


sdMI 5 

oVXIVXX . (J 


PHG 

X V,XX VJ 


489 


T/P515-3 


998 


853 


766 


889 


-95 


< 

o 


V/OVXIVXVJ. VJ 


60 

VJVJ 


4Q1 

It/1 


d 91fl-D19 

V.i iilU ulf 


843 


fl 708 


D 589 


D 759 


-149 


9 

i J 


sHMI 5 

oVXIVX X . U 


60 

VJVJ 




n ni8-D5i 

V.i UIO UtJl 


1 094 


n 996 


D 854 


D 999 


-1 59 


e 

: O 


psHK7 

CoVXXV i 


6fl 

VJVJ 


536 

UOVJ 


n 198-D34 

V.I XZjO UOt: 


883 


784 


D 659 


D 896 


-57 


8 

: O 


oVXIVXVJ. (J 


60 

VJVJ 


537 


O 098-043 

V.I VJZjO VJtiO 


1 DD9 

X . VJVJZj 


1 nm 

X . VJ XVJ 


D 978 


D 989 


-1 1 3 

X xo , 


n 

VJ 


K5 V 

XV(J V 


60 

VJVJ 


1 D88 

1UOO 


n 91 7-D55 

V.i Zv X 1 _ VJ<JtJ 


D 454 


n 78i 

U. 1 ox 


D 497 


D 74D 


1 7 

X f 




M3 5 V 

IVXO.O V 


60 

VJVJ 


1 1 64 


L 990-97 

XJ ZjZjVJ I 


D 654 


fl 868 

VJ .QUO 


D 988 

U. JOG) 


D 955 


-99 




K5 V 

XV O V 


1 00 

XVJVJ 


1 1 74 


LP 406-47 


924 

VJ . •JZj'X 


607 

VJ . VJVJ 1 


454 


659 


-111. 


7 


COVX1VXO. VJ 


200 a 


1481 


LP 711-32 

XJX 1 X X O — . 


554 

VJ . ' ) 'Jl 


708 

VJ . 1 VJO 


443 


681 

U. VJO X 


94 


. O 


sdM3 

oVXIVXO. VJ 


PHG 

X V,XX VJ 


1 749a 


TP 41 7-49 


D 689 

VJ. VJOi/ 


fl 498 

VJ .t:ZjO 


D 984 


D 499 


91 8 

Lid 


1 

. X 


CoVXlVXU . <J 


900 

^(VJVJ 


1 970 

X tJ 1 VJ 


LP 484-6 


878 


D 589 

VJ .UOi/ 


D 475 


D 657 


9 


s 


esrlM9 5 


200 


91 1 n 

_■ x xu 


LP 787-4 


D 369 


fl 747 


D 366 


D 653 


37 


. o 


M4 n V 

lVX i ±.VJ V 


1 00 

XVJVJ 


2497 


Gl 455 

V.I 1 t: • J ' J 


542 

VJ . Ot:Zj 


656 


400 


627 


36 


2 


srlM3 5 

oVXlVXO. iJ 


PHG 

X V,XX VJ 


2715 

Zj 1 X (J 


Gl 506 1 

V.I 1 OVJVJ . X 


994 

VJ . J 


1 019 


971 


975 


36 


n 

. VJ 


OVXX V 


60 

VJVJ 


9859 


TP 856-36 

XJX ouu ou 


n 651 


fl 736 


D 5D7 


D 795 


1 3 

XO, 


< 


DVllVXZ; . VJ 


XVXXVJ 


3061 

OVJVJ X 


T,P 5D2-32 

XJX <J VJ j-j tJ — 


D 737 

VJ . ( O 1 


D 497 

VJ . T.J 1 


0.292 


D 497 


65 

VJO . 


4 


COVXIVXO • VJ 


200 


3073 

OVJ ( o 


G 1 37-DD8 

V.I XO 1 IJVJO 


f) 989 

VJ . Ju J 


n 91 8 


D 869 


0.922 


-254. 




oVXXV i 


60 

VJVJ 


3084 


G D1 5-D26 

V.I UlU VJZj VJ 


f) 84D 


D 878 


D 724 


D 857 


-43 


« 




60 


o -i on 

31o9 


Lr ZZo-ZZ 


U.341 


U.63 f 


n on k 
U.3U5 


n £^7n 
U.o r 9 


-1U6. 


.1 


salvlo.U 


onn 
2UU 


3192 


Gl 871 


0.454 


0.768 


0.433 


0.706 


0. 


.6 


M3.5 V 


60 


3193 


G 169-007 


0.688 


0.878 


0.635 


0.829 


20. 


.1 


M1.0 V 


60 


3259 


LP 686-36 


0.995 


0.811 


0.657 


0.803 


-217, 





esdM0.5 


200 


3382 


LP 24-219 


0.937 


0.635 


0.485 


0.667 


-138. 


9 


esdM2.5 


200 a 


3409 


LP 141-1 


0.443 


0.645 


0.355 


0.587 


-80, 


3 


sdM4.5 


RHG 


3480 


LP 869-24 


0.458 


0.695 


0.363 


0.632 


12 


.9 


sdM4.0 


200 


3481 


LP 753-21 


0.782 


0.669 


0.467 


0.694 


75. 


.4 


sdM2.5 


200 


3548 


LP 695-96 


0.886 


0.839 


0.435 


0.660 


80. 


.6 


esdM3.0 


200 a 


3628 


LP 757-13 


0.908 


0.643 


0.576 


0.742 


10. 


.1 


esdM1.5 


100 
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Table 2 — Continued 



LHS Name Ti05 CaHl CaH2 CaH3 v rad class Source 



3684 


LP 518-12 


0.377 


0.733 


0.344 


0.637 


-75.3 


M4.5 V 


200 


3867 


G 067-030 


1.024 


0.997 


0.962 


0.982 


-108.8 


sdK 


60 


3957 


G 190-026 


0.727 


0.846 


0.648 


0.842 


-85.0 


M0.5 V 


60 


6304 


G 169-009 


0.587 


0.785 


0.505 


0.751 


-9.1 


M2.5 V 


60 




G 017-006 


0.781 


0.920 


0.702 


0.875 


8.9 


MO.O V 


60 




Sm 183 


0.880 


0.947 


0.867 


0.920 


-12.3 


gM 


100 




TVLM 832-42500 


0.186 


0.842 


0.171 


0.452 


72.1 


>M6.0 V 


200 



a Observed with August 1995 setup 



Table 3. Model Best Fits 



Star 


Class 


T <=// 


[m/H] 


Star 


Class 


T eff 


[m/H] 


Star 


Class 


T eff 


[m/H] 


LHS 410 


M1.0 V 


3700 


-1.0 


LHS 


307 


sdM0.5 


3700 


-1.5 


LHS 3259 


esdM0.5 


3700 


-2.0 


LHS 301 


M3.5 V 


3700 


-0.5 


LHS 


29 


sdMl.O 


3700 


-1.5 


LHS 192 


esdMl.O 


3600 


-2.5 


LHS 218 


M3.5 V 


3600 


-1.0 


LHS 


479 


sdMl.O 


3700 


-1.0 


LHS 364 


esdM1.5 


3600 


-2.5 


LHS 3192 


M3.5 V 


3600 


-0.5 


LHS 


178 


sdM1.5 


3600 


-1.0 


LHS 3628 


esdM1.5 


3600 


-1.5 


LHS 376 


M3.5 V 


3600 


-0.5 


LHS 


482 


sdM1.5 


3600 


-1.0 


LHS 161 


esdM2.0 


3600 


-2.0 


LHS 6304 


M2.5 V 


3600 


-1.0 


LHS 


491 


sdM1.5 


3600 


-1.5 


LHS 3382 


esdM2.5 


3500 


-1.5 


LHS 1088 


M3.5 V 


3500 


0.0 


LHS 


64 


sdM1.5 


3600 


-1.0 


LHS 1970 


esdM2.5 


3500 


-1.5 


LHS 2110 


M4.0 V 


3500 


0.0 


LHS 


320 


sdM2.0 


3600 


-1.0 


LHS 1174 


esdM3.0 


3500 


-1.5 


LHS 55 


M3.5 V 


3500 


0.0 


LHS 


216 


sdM2.0 


3600 


-1.0 


LHS 3548 


esdM3.0 


3500 


-1.5 


LHS 57 


M4.0 V 


3500 


-0.5 


LHS 


425 


sdM2.0 


3600 


-1.0 


LHS 453 


esdM3.5 


3400 


-2.5 


LHS 3684 


M4.5 V 


3500 


-0.5 


LHS 


2852 


sdM2.0 


3600 


-1.0 


LHS 375 


esdM4.0 


3400 


-2.0 


LHS 460 


M5.0 V 


3300 


0.0 


LHS 


20 


sdM2.5 


3600 


-1.0 


LHS 205a 


esdM4.5 


3400 


-2.0 


LHS 254 


M6.5 V 


3200 


0.0 


LHS 


3481 


sdM2.5 


3600 


-1.5 


LHS 3061 


esdM5.0 


3300 


-2.0 


TVLM 832-42500 


>M6.0 V 


3000 


0.0 


LHS 


272 


sdM3.0 


3600 


-1.5 


LHS 1742a 


esdM5.5 


3300 


-2.0 










LHS 


156 


sdM3.0 


3600 


-1.5 


















LHS 


1481 


sdM3.0 


3600 


-1.0 


















LHS 


2497 


sdM3.5 


3500 


-1.0 


















LHS 


3480 


sdM4.0 


3500 


-1.0 


















LHS 


3409 


sdM4.5 


3400 


-1.0 


















LHS 


3189 


sdM5.0 


3400 


-1.0 


















LHS 


407 


sdM5.0 


3400 


-1.5 


















LHS 


377 


sdM7.0 


3200 


-1.5 
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Table 4. Photometry 



LHS V B-V V-R R-I V-I vr (") a T (") A LK a My Source b 



12 


12.24 


1 


.45 


0. 


.90 





.90 


1. 


,80 


0.0359 


0.0031 


-0.08 


10.02 


Y,L 


20 


13.01 


1 


.68 


1, 


.03 


1 


.17 


2, 


.20 


0.0648 


0.0025 


-0.01 


12.07 


Y,L 


29 


8.84 


1 


.56 


0. 


.96 


1 


.00 


1. 


,96 


0.2583 


0.0065 


-0.01 


10.90 


Y,L 


35 


12.83 


1 


.72 


1, 


.25 


1 


.67 


2, 


.92 


0.1480 


0.0022 


-0.00 


13.68 


Y,L 


42 


12.23 


1 


.45 


1, 


.06 





.92 


1. 


,98 


0.0327 


0.0025 


-0.06 


9.80 


Y,E79 


44 


6.45 





.76 














0.1127 


0.0014 


-0.00 


6.71 


Y,S 


52 


9.43 





.84 


0. 


.52 





,49 


1. 


.01 










R 


53 


9.04 





.78 


0. 


.45 





.46 





.91 










R 


55 


12.90 


1 


.54 














0.0550 


0.0029 


-0.03 


11.60 


HD,S 


64 


13.30 


1 


.55 














0.0419 


0.0022 


-0.03 


11.41 


Y,S 


103 


14.16 






1. 


.15 


1. 


.47 


2 


.62 


0.0536 


0.0021 


-0.02 


12.81 


Y,WU 


104 


13.78 


1 


.34 





.81 


0. 


.91 


1. 


.72 










E79 


137 


13.36 


1 


.60 














0.0463 


0.0027 


-0.03 


11.69 


Y,S 


139 


15.05 


1 


.86 


1, 


.24 


1. 


.61 


2 


.85 


0.0485 


0.0044 


-0.09 


13.48 


Y,E79 


156 


14.89 


1 


.70 


0. 


.92 


1 


,25 


2 


,17 


0.0277 


0.0035 


-0.18 


12.1 


Y,E79 


161 


14.75 


1 


.55 


1, 


.01 





.96 


1. 


,98 


0.0260 


0.0048 


-0.42 


11.82 


Y,E79 


163 


13.06 


1 


.57 


1 


.12 


1 


.43 


2, 


,55 


0.0496 


0.0043 


-0.08 


11.54 


Y,B 


169 


14.13 


1 


.45 





.91 





.76 


1. 


,72 


0.0309 


0.0023 


-0.06 


11.58 


Y,L,E79 


170 


10.69 


1 


.22 


0. 


.80 





.70 


1. 


.50 


0.0348 


0.0025 


-0.05 


8.40 


Y,E79 


173 


11.12 


1 


.31 














0.0353 


0.0031 


-0.08 


8.86 


Y,S 


174 


12.75 


1 


.52 


0. 


.94 





.89 


1. 


.83 


0.0204 


0.0045 


-0.65 


9.30 


Y,DF89 


175 


9.91 





.65 














0.0218 


0.0022 


-0.11 


6.6 


Y,R 


178 


12.90 


1 


.55 


0. 


.95 


1 


.11 


2, 


.06 


0.0451 


0.0066 


-0.24 


11.17 


Y,E79 


182 


13.42 


1 


.57 














0.0234 


0.0025 


-0.12 


10.27 


Y,S 


185 


15.30 


1 


.79 


0. 


.98 





,85 


1 


.83 


0.0167 


0.0046 


L 


11.41 


Y,DF89 


192 


17.33 














1 


.98 


0.0102 


0.0008 


-0.06 


12.37 


M92 


205a 


18.93 














2 


.35 


0.0104 


0.0013 


-0.17 


14.02 


M92 


211 


14.11 


1 


.45 


0. 


.91 


1 


.05 


1. 


.96 


0.0188 


0.0029 


-0.28 


10.48 


Y,L 


216 


14.66 


1 


.62 


0. 


.99 


1 


.09 


2 


.08 


0.0306 


0.0030 


-0.10 


12.09 


Y,L 


218 


14.84 


1 


.57 


1 


.12 


1 


.33 


2. 


.45 


0.0296 


0.0030 


-0.11 


12.20 


Y,DF89 


232 


13.72 


1 


.15 














0.0148 


0.0025 


-0.34 


9.57 


Y,S 


236 


13.10 


1 


.33 














0.0187 


0.0018 


-0.10 


9.46 


Y,S 


241 


8.32 





.62 






















S 


254 


17.41 


1 


.75 


1, 


.90 


2 


.07 


3 


,97 










L 


272 


13.46 






1, 


.34 


1 


.36 


2, 


,70 










DF89 


276 


11.86 


1 


.43 


0. 


.93 


1 


.04 


1. 


.97 


0.0347 


0.0025 


-0.05 


9.56 


Y,B 
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Table 4 — Continued 



LHS 


V 


B-V 


V-R 


R-I 


V-I 


7T (") 


°* ( ) 




My 


Source 


3m 

GUI 


1 3 ^7 


1 ^7 


1 


08 


1 


39 


9 47 


n n c i4n 


n nn36 


n ni 


1 9 93 


V W88 F,7Q 

I . VV OOj-Cj / V 


3n7 

OU i 


1 5 99 

_1 ■ ) . Z, _■ 


1 56 










1 99 

1 . Cf Zi 


n ni 83 


n nn97 

U.UUZj i 


-n 95 

U. ZjO 


1 1 53 

1 1 . OO 


Y J,* 


320 


14.00 


1.54 





88 


1 


26 


2.14 


n n26n 


0.0034 


-9.19 


11.07 


Y,E79 


343 


13.87 


1.30 





87 





69 


1.56 


n.ni67 


0.0031 


-9.43 


9.98 


HD,L 


364 


14 61 


1 71 

1.11 


1 


03 





92 


1 95 


n n374 


n nn37 


-n in 

U. 1U 


1 9 47 

1 Zj .4 t 


Y DF99 


375 

O i U 


1 5 68 

J. o .uo 


1 87 


1 


08 


1 


12 


9 9n 

_■ . ZvU 


n n395 

U. UOi/O 


n nnm 

U.UU1U 


-n ni 

U. U 1 


1 3 66 

lO .UU 


T> 4 

1 lil 


376 

O i U 


1 5 nn 


1 64 

1 .U4 


1 


08 


1 


33 


9 41 
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Table 4 — Continued 



LHS 
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B-V 


V-R 


R-I 


V-I 


tt(") 


(") 




M v 


Source b 


3073 


13.69 


1.41 


0.84 


0.74 


1.58 


0.0187 


0.0034 


-0.41 


10.05 


Y,DF92 
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13.43 


1.43 


0.93 


0.90 
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0.0190 


0.0025 


-0.19 


9.82 
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0.0030 
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14.3 
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1.70 
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3192 
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1.11 
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0.0023 
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15.16 


1.85 
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2.82 


0.0500 


0.0013 
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0.0008 
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13.48 
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0.0007 
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0.0006 
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12.12 
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0.0125 
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RA 
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0.0008 
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Y,S 


4037 


13.52 
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0.47 
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6304 


14.07 




1.05 


1.25 


2.30 










W84 



a L indicates the Lutz-Kelker correction is undefined but < —0.80. 

b References for the trigonometric parallaxes and photometry are: B - Bessell (1990), 
CNS3 - Gliese and Jahreiss (1991), DF89 - Dawson and Forbes (1989), DF92 - Dawson 
and Forbes (1992), E79 - Eggen (1979), E87 - Eggen(1987), HD - Harrington and Dahn 
(1980), L - Leggett (1992), M92 - Monet et al. (1992), R - Ryan (1992), RA - Ruiz and 
Anguita (1993), S - SIMBAD, T - Tinney et al. (1995), W84 - Weis (1984), W88 - Weis 
(1987), W91 - Weis (1991), WU - Weis and Upgren (1982), Y - Van Altena et al. (1991) 

°Includes 0.75 magnitude correction for unresolved companion 
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Table 5. Space Velocities 
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Table 6. Velocity Statistics 
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sdM 
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5 
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61 


strong CaHl 


24 


-38 


61 


44 


30 


weak CaHl 
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48 


30 


24 


Disk (RHG) 
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43 
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I. HBO 



Fig. 1. — The CaH indices as a function of Ti05. Stars classified as sdM appear as filled squares, 
esdM as filled triangles, sdK as filled circles, and spectroscopic non-subdwarfs as open triangles. 
In the CaHl diagram we also plot the single stars within eight parsecs as open squares. All other 
figures use the same symbols. Representative error bars (±0.03) appear in the lower right. 
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Fig. 4.— The CaH3/CaH2 ratio as function of CaH2. The sdM and esdM lie below the RHG 
Population I relation. The RHG stars appear as dots. 




Fig. 5. — A model atmosphere fit to the Population I star Gl 402, an M 4.0 V standard. Note that 
the overly strong atomic lines are an artifact of the computation technique. 
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Fig. 7. — The observed HR diagram. The spectroscopic extreme subdwarfs (esdM) are filled 
triangles, the intermediate subdwarfs (sdM) are filled squares, and the subdwarf candidates that 
were spectroscopically indistinguishable are open triangles. The open squares are the single stars 
within eight parsecs with good parallaxes - note the step at V-I ~ 2.8 which is discussed in HGR. 
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Fig. 8. — The HR diagram for B-V colors. The three pointed stars are nearby G and K dwarfs. 
Note that the blue K-subdwarfs lie below the disk sequence, but the sdM and esdM sequences cross 
the disk main sequence and are actually brighter at the reddest colors. 



_ 41 - 



_ 1 1 1 

— 




I 


1 


1 




1 


1 


- 
















— 








< 










□ 




□ 










- 


□ 












- 




„ □ 
□ 

□ 




□ 


■ 








- 




□ 














□ 
□ 














n 


















< 




■ 










— 


< 




■ 




■ 




— 


— 


□ 














~~ n 

□ n 
□ □ 


■ 

■ 


■ 

I 

A 








< 




□ □ 

n 

<■□ 

< 


■ 
■ 


■ << 












□ 


















■ 






4 










■ 
















i 


1 


1 


1 






1 



CD 



00 



CM 



CM 

c\i 



A-a 



Fig. 9. — V-I vs. B-V. The sdM are redder in B-V at a given V-I color than the disk stars. The 
esdM show a larger offset than the sdM 
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Fig. 10. — The observed HR diagram with theoretical sequences. Metallicities of [m/H] = 0.0 
(dotted line ), -0.5 (short dashed line), and -1.5 (long dashed line) from Baraffe et al. (1995). 
A more recent solar metallicity model (solid line) computed by Baraffe and Chabrier using the 
NextGen model atmospheres is also shown. The latter model does not match the observed disk 
stars, but the agreement of the previous models is quite good. 
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Fig. 11. — The observed HR diagram with the cluster main sequences fit by Santiago et al.(1996) 
using HST observations. From left to right, the clusters are M15 (solid line, [Fe/H] = —2.26), 
u Cen (dotted, -1.6), 47 Tuc (short dashed, -0.6), N2420 (long dashed, -0.45), and N 2477 (long 
dash-dotted, 0.00). Note that the position of the cluster sequences agrees with our spectroscopic 
metallicity estimates. The slopes of the cluster sequences do not disagree with the local subdwarfs 
in the region of color overlap when using our classifications. The discrepancy of the reddest stars 
in the metal rich clusters may be due to color terms, as well as the inadequacy of a linear fit for 
the main sequence. 



